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LONG-TERM  GOALS 

This  study  contributes  to  our  long-term  efforts  toward  understanding: 

•  Mixed  layer  dynamics 

•  Processes  that  communicate  atmospheric  forcing  to  the  ocean  interior. 

OBJECTIVES 

Existing  high  resolution  regional  models  typically  resolve  the  mean  vertical  structure  of  the  upper 
ocean  boundary  layer.  Physically-based  parameterizations  of  vertical  fluxes  make  it  possible  to 
account  for  subgrid  mixing  at  length  scales  smaller  than  the  layer  depth,  but  no  specialized 
parameterization  is  used  to  represent  the  dynamics  of  horizontal  mixing  below  the  0(1  )km  -  0(10 )km 
resolution  scale.  We  aim  to  determine  the  physical  limitations  of  subgrid  parameterization  on  these 
scales.  This  project  addresses  the  following  questions: 

•  What  physics  govern  horizontal  and  vertical  mixing  in  the  presence  of  horizontal  variability  on 
the  1-10  km  scale? 

•  What  is  the  relative  importance  of  horizontal  and  vertical  mixing  in  determining  the  structure 
of  the  boundary  layer? 

•  How  well  can  existing  parameterizations  simulate  vertical  and  horizontal  mixing? 

•  What  physics  should  be  included  to  improve  parameterizations? 

APPROACH 

An  adaptive  measurement  program  employed  acoustically-tracked,  neutrally  buoyant  Lagrangian 
floats  and  a  towed,  undulating  profiler  to  investigate  the  relative  importance  of  vertical  and  horizontal 
mixing  in  governing  boundary  layer  structure  in  the  presence  of  0(1  km)  scale  horizontal  variability. 
Remotely  sensed  sea  surface  temperature  and  ocean  color,  combined  with  rapid,  high-resolution  towed 
surveys  and  model  results  guide  float  deployments  to  key  locations  within  fronts.  Synoptic,  high- 
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resolution  surveys  followed  Lagrangian  float  drifts  to  characterize  three-dimensional  variability  within 
the  span  of  a  model  grid  points.  Acoustic  tracking  allowed  towed  surveys  to  follow  floats  and 
geolocated  all  observational  assets  for  later  analysis.  Measurements  characterized  boundary  layer 
turbulence  and  facilitate  detailed  separation  of  vertical  and  horizontal  processes. 

A  turbulence-resolving  Large  Eddy  Simulation  (LES)  was  used  to  model  the  dynamics  of  vertical  and 
horizontal  mixing  in  a  domain  volume  corresponding  to  a  regional  model’s  horizontal  gridscale  and  set 
in  the  translating  Lagrangian  reference  frame  of  the  float/survey  observations.  The  observations  will 
provide  realistic  initial  and  time-dependent  boundary  conditions  and,  in  particular,  time-dependent 
lateral  boundary  conditions  will  be  determined  from  rapid  surveys. 

Quantitative  one-to-one  statistical  comparisons  between  LES  results  and  the  float  and  survey 
observations  will  be  made.  This  product  will  have  direct  application  to  assessing  regional  model 
subgrid  parameterizations. 

WORK  COMPLETED 

The  2006  Measurement  Program:  Observations  in  the  California  Current 
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Figure  1.  Float  drift 
(black)  and  towed  profiler 
(magenta)  tracks  plotted 
over  remotely  sensed  sea 
surface  temperature.  The 
bathymetry  rises  from  the 
abyss  near  the  eastern 
margin  of  the  chart. 


The  first  of  two  cruises  associated  with  this  AESOP  effort  took  place  from  R/V  Roger  Revelle,  16  July 
-  8  August  2006  off  the  California  coast.  Operations  began  with  instrument  testing  and  acoustic 
tracking  refinement  in  the  Channel  Islands  and  Santa  Barbara  channel.  Following  this,  two  drifting 
surveys  focused  on  a  zonally  oriented  front  located  west  of  the  continental  rise  off  of  San  Luis  Obispo 
(Fig.  1).  A  third  drift  followed  the  southward  flow  associated  with  a  strong  meridionally  oriented  front 
(Fig.  2).  Sections  occupied  prior  to  float  deployment  exhibit  T-S  characteristics,  small  pycnostads  and 
optical  signatures  that  suggest  active  subduction  of  cold-side  waters  into  the  region  below  the  warm- 
side  mixed  layer  base  (Fig.  2).  Informed  by  float  behavior  inferred  from  model  results  analyzed  before 
the  cruise  and  by  the  high-resolution  towed  profiler  section  occupied  immediately  prior  to  float 
deployment,  we  selected  a  site  intended  to  maximize  the  probability  of  observing  subduction.  Our  first 
few  deployments,  during  periods  of  weak  wind,  found  little  evidence  for  subduction  driven  by 
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Figure  2.  (left)  Front-following  drift  with  line  colors  as  defined  in  Fig.  1.  (right)  Potential  density 
(contours)  and  660  nm  beam  attenuation  (colors)  from  the  section  used  to  choose  the  float  deployment 
site. 

submesoscale  processes  alone.  Our  final  two  deployments,  spanning  a  period  of  rising,  sustained  and 
then  relaxing  winds,  clearly  showed  subduction  and  restratification  as  the  front  relaxed  after  the  wind 
dropped. 

The  2007  Measurement  Program:  Observations  in  the  Kuroshio  Extension 

Building  on  lessons  learned  working  off  the  California  coast  in  2006,  the  2007  field  effort  focused  on 
the  strong  fronts  and  submesoscale  features  associated  with  the  Kuroshio  extension.  As  during  the 
previous  cruise,  near-realtime  remotely  sensed  images  and  numerical  results  guided  sampling  to 
promising  features.  This  resulted  in  three  surveys,  described  below. 

The  first  deployment  (Fig.  3)  was  placed  nearby  in  the  Kuroshio  extension  third  meander  and  focused 
on  a  meridionally-oriented  section  of  front  forced  by  southerly  (up-front)  winds.  We  anticipated,  and 
found,  that  the  upfront  winds  would  drive  frontalysis,  spreading  light,  warm-side  waters  over  denser 
cold-side  waters.  We  hoped  that  blocking  of  Ekman  transport  by  the  fronts’  relative  vorticity  field 
would  cause  downwelling.  However,  the  front  remained  weak  during  the  roughly  four  day  survey  with 
the  float  experiencing  no  downwelling  events. 

The  second  survey  (Fig.  4)  focused  on  a  steep  southward  meander  that  NLOM  forecasts  predicted 
would  detach  within  a  few  days  time.  Strong  northerly  winds  provided  frontogenic  (downwind) 
forcing  along  the  western  side  of  the  loop  and  frontolytic  (upwind)  forcing  over  the  eastern  side.  The 
pre-deployment  survey  revealed  strong  fronts  marking  both  sides  of  the  feature,  with  the  deepest 
mixed  layers  along  the  western  margin.  We  deployed  a  float  at  the  strongest  gradient  region  of  the 
front  on  the  western  side,  calculating  that  it  would  transit  over  the  ‘U’-shaped  path  over  3-4  days 
contrasting  the  response  to  down-front  and  up-front  wind  forcing.  Over  four  days  of  weakening  winds, 
the  fronts  marking  both  sides  of  the  meander  became  more  diffuse,  with  a  buoyant  cap  restratifying  the 
mixed  layer  across  the  entire  feature. 
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Figure  3.  First  Kuroshio  survey.  Colors  show 
MODIS  SST  image.  Inset  shows  the  survey. 
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Figure  4.  Second  Kuroshio  survey  in  a  detaching 
meander  -  “The  Sock” 


40°N 


38°N 

jse-N 


34°N 


32°N 


30°N 


138°E  141°E  144°E  147°E  150°E 

longitude 


Figure  5.  The  third  survey 
of  the  2007  cruise,  focusing 
on  the  first  section  of  the 
Kuroshio  extension.  The 
survey  (lower  insert) 
followed  the  front  as  it  is 
advected  westward  from  its 
creation  off  the  Japanese 
coast  by  the  convergence  of 
the  cold  Oyashio  and  warm 
Kuroshio  currents. 


The  final  survey  focused  on  the  strong,  nearly  linear ,  sharp  and  zonally-oriented  section  of  Kuroshio 
front  extending  offshore  from  Japan  (Fig.  5).  Forcing  by  northwesterly  (downfront)  winds  provided  an 
ideal  combination.  The  pre-deployment  section  revealed  intense  cyclonic  relative  vorticity  and  optical 
signals  indicative  of  recent  subduction  in  strong  gradient  region.  Over  a  7-day  period,  Triaxus 
repeatedly  occupied  cross-front  sections  following  the  drifting  float,  attempting  to  characterize  a 
region  marked  by  two  outcropping  isopycnals.  Results  from  this  survey  are  described  in  detail  below. 
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Data  Processing  and  Analysis 

Triaxus  towed  profiler,  float  and  shipboard  data  have  been  processed  for  both  cruises.  Triaxus- 
mounted  ADCP  data  have  been  processed  for  both  deployments.  Analysis  efforts  have  focused  on 
restratification  observed  at  a  front  in  the  California  Current  and  on  the  evolution  of  a  particularly 
sharp,  wind-forced  from  in  the  Kuroshio.  The  Kuroshio  analysis  resulted  in  a  paper  published  in 
Science,  while  the  California  Current  analysis  is  still  in  progress.  This  work  involves  collaboration 
with  Dr.  Andrey  Shcherbina  (APL-UW)  and  Dr.  Leif  Thomas  (Stanford). 

RESULTS 

Please  see  the  attached  paper,  published  in  Science. 

IMPACT/APPLICATION 

None. 

TRANSITIONS 

None. 

RELATED  PROJECTS 

SeaSoar  and  Doppler  Sonar  Spatial  Survey  of  Internal  Tide  Generation  and  Mixing,  Shaun  Johnston 
and  Daniel  Rudnick. 

Lagrangian  Studies  of  Lateral  Mixing  (a  component  of  the  Lateral  Mixing  DRI).  Craig  Lee  and  Eric 
D’Asaro. 
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Enhanced  Turbulence  and  Energy 
Dissipation  at  Ocean  Fronts 

Eric  D'Asaro/*  Craig  Lee/  Luc  Rainville/  Ramsey  Harcourt/  Leif  Thomas2 

The  ocean  surface  boundary  layer  mediates  air-sea  exchange.  In  the  classical  paradigm  and 
in  current  climate  models,  its  turbulence  is  driven  by  atmospheric  forcing.  Observations  at  a 
1-kilometer-wide  front  within  the  Kuroshio  Current  indicate  that  the  rate  of  energy  dissipation 
within  the  boundary  layer  is  enhanced  by  one  to  two  orders  of  magnitude,  suggesting  that  the 
front,  rather  than  the  atmospheric  forcing,  supplied  the  energy  for  the  turbulence.  The  data 
quantitatively  support  the  hypothesis  that  winds  aligned  with  the  frontal  velocity  catalyzed  a 
release  of  energy  from  the  front  to  the  turbulence.  The  resulting  boundary  layer  is  stratified  in 
contrast  to  the  classically  well-mixed  layer.  These  effects  will  be  strongest  at  the  intense  fronts 
found  in  the  Kuroshio  Current,  the  Gulf  Stream,  and  the  Antarctic  Circumpolar  Current,  all  of 
which  are  key  players  in  the  climate  system. 


Although  the  basic  characteristics  of  ocean 
circulation  have  been  well  known  for 
many  decades,  a  detailed  understanding 
of  its  energetics  has  emerged  only  recently  (7). 


The  energy  sources  are  well  understood:  Wind 
stress  acting  on  surface  currents  (or  “wind-work”), 
particularly  in  the  Southern  Ocean,  is  the  dom¬ 
inant  energy  source,  with  little  net  input  from 


heating/cooling  or  precipitation/evaporation.  The 
energy  sinks,  however,  are  less  well  understood. 
Energy  dissipation  requires  a  cascade  of  energy 
through  nine  orders  of  magnitude,  from  the  size 
of  the  ocean  to  the  centimeter  scales  of  viscous 
dissipation.  A  cascade  of  processes  supports  this 
flux.  Instabilities  of  the  large-scale  circulation 
lead  to  the  generation  of  a  rich  field  of  eddies 
with  typical  scales  of  100  km  at  mid-latitudes. 
The  dynamics  of  these  eddies  are  highly  con¬ 
strained  by  Earth’s  rotation  such  that  their  cur¬ 
rents  are  nearly  geostrophic  (that  is,  the  flow  is 
governed  by  a  balance  between  Coriolis  and  hor¬ 
izontal  pressure  forces).  A  turbulent,  geostrophic 
eddy  field  tends  to  flux  energy  to  larger  rather 
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Fig.  1.  (A)  Sea  surface  temperature  (SST)  map 
(37)  of  experimental  area  with  sea  surface  height 
(SSH)  (8)  shown  as  contours  (0.1-m  interval).  The 
thick  black  line  denotes  the  ship  track.  Blue  and 
red  arrows  mark  sections  shown  below  in  blue  and 
red  boxes  (C  to  F),  respectively.  L  and  H  mark  low- 
and  high-pressure  centers,  respectively.  (B)  Triaxus 
and  float  data  around  SF2  [blue  box  (C  and  E)], 
with  Triaxus  density  shown  as  a  colored  curtain 
hanging  beneath  the  ship  track  and  float  trajectory 
shown  as  a  ribbon  colored  by  density  beneath  the 
surface  track  (dashed  line).  Sections  show  along- 
front  velocity  [(C)  and  (D)]  and  total  PV  [(E)  and  (F)] 
for  SF2  (left,  blue  box)  and  for  a  later  section  (right, 
red  box).  Black  contours  plot  potential  density  in 
each  section,  with  the  thick  line  marking  o0  =  25.0 
kg  m-3  and  thinner  lines  plotted  at  0.05  intervals. 
Gray  lines  denote  the  Triaxus  profile  track. 


A 


0 

73 

=5 

4— > 

CO 


B 


140 


o 

10 

20 

30 

40 

50 


142 


144  146 

longitude 


East  [km]  10 


148 


23°C 

22 

21 

20 

19 

18 

17 

16 

15 

14 


125.0 


J24.8 


124.6  £• 


(D 


jo  9  82. 

</)  O 

do 

'-5  J? 


318 


15  APRIL  2011  VOL  332  SCIENCE  www.sciencemag.org 


Downloaded  from  www.sciencemag.org  on  March  19,  2012 


than  smaller  scales,  thus  providing  no  obvious 
path  to  dissipation.  Recent  simulations  (2-4)  with 
very-high-resolution  models  suggest  a  new  path 
from  the  eddy  field  toward  dissipation  through  the 
formation  of  submesoscale  fronts,  regions  of 
strong  lateral  gradient  in  the  upper  ocean,  with 
horizontal  scales  of  1  to  10  km.  Instabilities  of 
these  fronts  could  then  cascade  energy  from  the 
frontal  scale  to  dissipation. 

The  simulations  also  suggest  that  the  surface 
“mixed-layer”  of  the  ocean  is  greatly  modified  in 
the  presence  of  fronts.  For  one,  the  boundary 
layer  is  stratified,  not  mixed,  and  deepens  by  the 
action  of  turbulent  motions  that  derive  part  of 
their  energy  from  the  frontal  circulation  as  op¬ 
posed  to  atmospheric  forcing  (4, 5).  This  is  a  shift 
from  the  classical  paradigm  of  a  surface  bound¬ 
ary  layer  driven  by  the  atmosphere,  with  implica¬ 
tions  for  climate  dynamics.  The  surface  boundary 
layer  is  the  mediator  for  air-sea  interaction  and 
therefore  influences  processes  that  play  an  inte¬ 
gral  role  in  the  climate  system  such  as  the  oceanic 
sequestration  of  carbon  and  the  subduction,  or 
transfer,  of  heat,  salt,  and  dissolved  gasses  from 
the  ocean’s  surface  to  its  interior.  Oceanic  se¬ 
questration  of  carbon  and  subduction  occur  to  a 


large  degree  in  the  proximity  of  the  ocean’s  main 
currents:  the  Gulf  Stream,  the  Antarctic  Circum¬ 
polar  Current,  and  the  Kuroshio  Current  (6,  7). 
These  currents  are  regions  of  strong  lateral  den¬ 
sity  gradients;  therefore,  their  surface  boundary 
layers  could  be  substantially  affected  by  frontal 
dynamics.  Here,  we  present  experimental  evidence 
from  a  front  in  the  Kuroshio  showing  that  this  is 
indeed  the  case. 

Evolution  of  the  front.  Measurements  were 
taken  from  18  to  21  May  2007  (days  137  to  140) 
near  the  start  of  the  Kuroshio  extension  off  the 
coast  of  Japan.  Here,  the  cold  subpolar  gyre  wa¬ 
ters  of  the  Oyashio  Current  meet  the  warm  sub¬ 
tropical  gyre  waters  of  the  Kuroshio  Current  to 
form  the  Kuroshio  front  (Fig.  1  A).  The  region  is 
rich  in  eddies,  as  illustrated  by  the  sea  surface 
height  contours  in  Fig.  1A  (5).  Measurements 
focused  on  an  exceptionally  sharp  front  [the 
“sharpest  front”  (SF),  blue  arrow  in  Fig.  1]  formed 
by  a  strongly  confluent  flow  between  two  eddies 
acting  on  the  contrast  between  the  warm  Kuroshio 
water  and  the  cold  Oyashio  water.  Measurements 
were  made  by  deploying  a  subsurface,  neutrally 
buoyant  Lagrangian  (water-following)  float  (9) 
in  the  frontal  region,  acoustically  tracking  it  from 
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Fig.  2.  (A)  SST  anomaly  (relative  to  SST  along  the  front)  from  ship  surveys  plotted  relative  to  float 
position.  Vertical  axis  is  cross-frontal  distance.  (B)  Depth-time  section  of  along-front  velocity  following  the 
float  with  potential  density  contoured  at  0.2-kg  m-3  intervals  where  available  in  the  upper  150  m.  (C) 
Wind  stress  (22,  23)  magnitude  and  along-front  component  in  pascals.  (D)  Vertical  shear  (Uz)  of  zonal 
velocity  as  a  function  of  depth  and  time  along  the  front. 
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the  ship,  and  surveying  a  ~5-km  region  around 
the  float  using  the  Triaxus  profiling  vehicle  towed 
behind  the  R/V  Melville  (10).  An  example  of  the 
sampling  strategy  is  shown  in  Fig.  IB.  The  float 
provided  a  reference  frame  moving  with  the  fron¬ 
tal  water,  so  that  changes  in  frontal  structure  can 
be  interpreted  as  attributable  to  its  temporal  evo¬ 
lution.  The  vertical  motion  of  the  float  within 
the  boundary  layer  also  provided  estimates  of 
the  turbulence  intensity  and  dissipation  rate.  Mea¬ 
surements  included  temperature,  salinity,  and 
pressure  on  both  platforms  and  velocity  profiles 
obtained  by  combining  ship  and  Triaxus  data 
(fig.  S3)  (10). 

Confluent  flow  (evident  in  Fig.  1A)  concen¬ 
trates  the  large-scale  north-south  temperature  and 
salinity  gradients  into  a  smaller  region,  thus  form¬ 
ing  the  SF.  During  day  137,  estimates  from  the 
large-scale  velocity  field  derived  from  satellite 
altimetry  show  a  north-south  convergence  of  north 
velocity  dv/dy  <  0  (v,  northward  velocity;  y,  north¬ 
ward  distance)  and  an  east-west  divergence  of  east 
velocity  du/dx  >  0  (u,  eastward  velocity;  x,  east¬ 
ward  distance),  both  of  magnitude  ~0.5  x  105  s-1. 
Estimates  from  direct  velocity  measurements  (fig. 
S3)  (10)  along  the  front  and  ±3  km  to  either  side 
confirm  the  north-south  convergence  and  the  ac¬ 
celeration  of  the  along-front  velocity  (du/dx  ~ 
-dv/dy  - 1 .2  ±  0.7  x  10-5  s-1).  The  observed  front 
thinned  laterally  by  about  a  factor  of  2  during  this 
period  (Fig.  2 A),  a  rate  close  to  that  predicted 
from  these  velocity  gradients,  thus  forming  and 
maintaining  the  SF  through  day  137.  Both  con¬ 
fluence  components  dropped  to  zero  by  day  138.5 
and  then  become  divergent  with  a  magnitude  of 
0.5  ±  0.5  x  10-5  s'1.  Thus,  the  SF  is  a  transient 
region  of  strong  density  gradient  generated  by 
strong  local  confluence  and  embedded  within  the 
larger-scale  Kuroshio  front. 

A  section  across  the  SF  (Fig.  1C)  shows  it  to 
be  less  than  1  km  wide  in  surface  density  and 
-20  m  deep  (11).  A  surface-velocity  maximum, 
the  frontal  jet,  is  found  on  the  warm  side  of  the 
front  (Fig.  1C).  The  shear  below  this  jet  extends 
across  the  frontal  region  and  has  a  substantial 
geostrophic  component,  but  with  a  shear  mag¬ 
nitude  that  is  roughly  half  that  expected  geo- 
strophically  from  the  horizontal  density  gradient. 
A  similar  section  taken  2  days  and  220  km  down¬ 
stream  (Fig.  ID)  shows  a  much  wider  (4  km)  and 
deeper  (60  m)  frontal  zone.  The  large  horizontal 
density  gradients  across  the  front  have  eroded, 
although  the  net  contrast  remains  about  the  same, 
and  the  volume  of  water  with  intermediate  prop¬ 
erties  has  increased  (Fig.  1,  C  and  D).  Thus,  the 
decay  of  the  SF  appears  to  be  due  to  a  combi¬ 
nation  of  local  difluence  and  mixing. 

Turbulence  and  mixing  at  the  front.  Mea¬ 
surements  collected  by  the  Lagrangian  float  at 
the  SF  quantify  the  rate  of  turbulent  mixing.  The 
float  interspersed  periods  of  Lagrangian  drift  with 
profiles  and  surfacing  for  communication.  During 
the  16  drifts,  the  float  was  water-following  in 
three  dimensions,  with  its  vertical  velocity  mea¬ 
suring  the  vertical  velocity  of  the  water.  During 
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profiles  and  surfacings,  the  float  purposefully 
moved  relative  to  the  water  and  was  thus  not  La- 
grangian.  The  float  sampled  the  SF  during  two 
of  the  drifts,  which  will  hereafter  be  called  SF1 
and  SF2.  During  these  drifts,  especially  SF2,  the 
float  was  exactly  at  the  front,  as  indicated  by:  (i) 
its  track  (Fig.  IB),  which  follows  the  frontal  in¬ 
terface  mapped  by  the  Triaxus;  (ii)  its  density 
(Fig.  3B),  which  is  intermediate  between  that  of 
the  warm  and  cold  sides  of  the  front;  and  (iii)  its 
location  at  the  maximum  in  horizontal  density 
gradient  (Figs.  IB  and  3 A).  During  these  drifts, 
the  float  repeatedly  cycles  across  the  upper-ocean 
boundary  layer  (Fig.  3E),  tracing  the  trajecto¬ 
ries  of  boundary-layer  water  parcels  and  thus 
measuring  their  vertical  velocity.  Numerous  other 
measurements  in  both  convection  (72,  13)  and 
wind-forced  boundary  layers  ( 14-1 7)  confirm  that 
the  vertical  motion  of  these  floats  is  due  to  upper- 
ocean  turbulence. 


What  is  the  energy  source  for  this  turbulence? 
Numerous  float  observations  of  the  average  ver¬ 
tical  kinetic  energy  (VKE)  in  the  upper-ocean 
boundary  layer  show  a  remarkably  good  correla¬ 
tion  between  VKE  and  the  10-m  wind  speed 
(Fig.  3C).  This  correlation  indicates  that  wind  is 
usually  the  major  source  of  energy  for  upper- 
ocean  turbulence  (18).  However,  the  VKE  during 
SF2  (gray  in  Fig.  3C)  lies  far  (a  factor  6)  above 
these  values,  suggesting  that  wind  cannot  explain 
the  anomalously  high  turbulence  levels  at  SF2. 

In  Fig.  3D,  the  anomalously  high  turbulence 
at  SF2  is  shown  in  terms  of  energy  dissipation 
rate;  that  is,  the  flux  of  energy  through  the  tur¬ 
bulence.  The  average  dissipation  in  the  boundary 
layer  is  estimated  from  the  frequency  spectra  of 
float  vertical  acceleration  (Fig.  3D,  black  cir¬ 
cles,  and  fig.  S4)  ( 10 , 19)  and  plotted  as  depth- 
integrated  dissipation  by  multiplying  by  the 
boundary-layer  depth  (Fig.  3E)  (20).  Dissipation 


rate  and  energy  are  closely  linked  so  that  a  sec¬ 
ond  estimate  with  higher  time  resolution  can  be 
formed  from  the  VKE  (Fig.  3D,  solid  and  dashed 
gray  traces  and  open  circles)  (21).  The  average 
dissipation  at  SF2  rises  by  a  factor  of  ~10  to  20 
above  the  nonfrontal  values  with  the  VKE  esti¬ 
mate  suggesting  that  even  higher  values  occurred 
just  before  the  measurements  began. 

Cooling  of  the  ocean  by  the  atmosphere  drives 
boundary-layer  turbulence  with  an  average  dis¬ 
sipation  rate  given  by  the  surface  buoyancy  flux 
(12).  At  SF2,  very  weak  cooling  occurs  (22,  23), 
with  a  buoyancy  flux  of  less  than  0.3%  of  the 
measured  dissipation.  Atmospheric  cooling  can¬ 
not  explain  the  anomalously  high  turbulence  lev¬ 
els  at  SF2. 

In  summary,  SF2  is  unremarkable  in  wind 
stress  (Fig.  2C)  or  in  velocity  (Fig.  2B),  but  is 
highly  anomalous  in  turbulence  level  (Fig.  3,  C 
and  D)  and  in  lateral  gradient  (Fig.  2A). 


Fig.  3.  (A)  Horizontal  buoyancy  gradient,  propor¬ 
tional  to  minus  the  density  gradient  (represented 
by  color),  computed  from  R/V  Melville's  hull- 
mounted  temperature  sensor  plotted  as  a  function 
of  time  and  potential  density.  The  black  line 
denotes  the  float  trajectory,  whereas  the  white  line 
plots  the  position  of  maximum  buoyancy  gradient. 
(B)  Density  profiles  at  SF2  from  float  (blue)  and 
Triaxus  (gray)  used  in  Fig.  IB.  Two  Triaxus  profiles 
are  highlighted  (black).  (C)  Vertical  velocity 
variance  from  floats  (14-17)  for  wind-forced 
upper-ocean  boundary  layers  as  a  function  of  wind 
speed.  Confidence  intervals  (95%)  were  computed 
as  in  (14).  SF2  lies  well  above  all  other  data.  U10/ 
the  wind  speed  at  10-m  height.  (D)  Time  series  of 
boundary-layer-integrated  dissipation  estimated 
from  the  float-acceleration  spectrum  (black  solid 
circle),  float  VKE  averaged  over  each  float  drift 
(open  circles),  1-hour-averaged  float  VKE  (gray 
line  with  dashed  line  interpolating  between  drifts), 
EBF  computed  at  float  averaged  over  drifts  (red 
circles),  and  EBF  computed  at  maximum  density 
gradient  (red  line,  1.5-hours  running  average).  (E) 
Float  depth  during  Lagrangian  drifts  (yellow  filled 
areas)  and  boundary-layer  depth  (black  bars) 
estimated  as  twice  the  mean  float  depth. 
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Potential  vorticity  and  frontal  instability. 

We  hypothesize  that  a  flux  of  energy  from  the 
front  itself  accounts  for  the  enhanced  turbulence 
levels  at  SF2.  The  boundary  layer  at  SF2  is  stably 
stratified  (Fig.  3B)  yet  highly  sheared  in  the  ver¬ 
tical  direction  due  to  the  presence  of  a  strong  jet 
along  the  front  (Fig.  2,  B  and  C).  This  latter  con¬ 
dition  makes  the  flow  potentially  susceptible  to 
symmetric  instability  (SI)  ( 4),  which  extracts  ki¬ 
netic  energy  from  the  geostrophic  frontal  jet.  The 
Ertel  potential  vorticity  (PV)  (24)  is  the  key  quan¬ 
tity  for  diagnosing  this  instability;  a  flow  is  un¬ 
stable  to  SI  when  the  PV  is  negative  (25).  PV  can 
become  negative  due  to  the  combination  of  a  suf¬ 
ficiently  strong  vertical  shear  and  lateral  density 
gradient  and  a  sufficiently  weak  vertical  density 
gradient.  These  conditions  can  occur  within  the 
boundary  layer  of  a  strong  front,  with  the  front 
providing  the  shear  and  lateral  gradient  and  the 
boundary  layer  having  a  reduced  stratification. 
Simulations  ( 4 ,  5)  indicate  that  under  these  con¬ 
ditions  SI  will  grow,  become  unstable  to  second¬ 
ary,  smaller-scale  instabilities  (26),  and  feed  a 
turbulent  cascade  to  dissipation,  resulting  in  a 
fully  turbulent  boundary  layer  drawing  its  energy 
from  the  front. 

We  used  velocity  and  density  data  taken  by 
the  ship  to  evaluate  the  PV  on  each  of  the  nearly 
100  crossings  of  the  front  (fig.  S5)  (JO,  24).  We 
found  negative  PV  near  the  surface  (Fig.  IE)  for 
0.2  days  at  SF2  and  nowhere  else  (Fig.  IF).  The 
front  at  SF2  is  therefore  unstable  to  SI,  suggest¬ 
ing  that  the  turbulence  at  SF2  is  drawing  energy 
from  the  frontal  shear. 

The  simulations  indicate  that  SI  at  a  front 
occurs  when  the  wind  blows  perpendicular  to  the 
frontal  gradient  (27,  28),  which  is  typically  in  the 
direction  of  the  frontal  velocity  (Fig.  4).  Such  a 


“down-front”  wind  drives  a  net  transport  of  water 
perpendicular  to  the  frontal  jet  to  carry  heavy 
water  across  the  front,  from  the  cold  side  to  the 
warm  side.  This  Ekman  transport  advects  heavy 
water  over  light  water,  reducing  the  stratification, 
and  thus  reducing  the  PV  and  promoting  SI.  The 
Ekman  buoyancy  flux  (EBF)  (27),  computed  from 
the  product  of  the  down-front  wind  stress  (Fig. 
2C)  and  the  cross-frontal  density  gradient  (Fig.  3  A), 
is  a  measure  of  this  effect.  Simulations  (4)  suggest 
that  turbulence  in  a  frilly  developed  boundary  layer 
of  depth  H  and  driven  by  down-front  winds  extracts 
kinetic  energy  from  the  frontal  jet  at  a  depth- 
integrated  rate  given  by  H(EBF)1 2  and  dissipates 
it  within  the  boundary  layer.  This  quantity  (29) 
peaks  at  SF2  (Fig.  3D,  red)  with  a  value  compara¬ 
ble  to  the  measured  dissipation  rate,  thus  provid¬ 
ing  quantitative  evidence  supporting  the  hypothesis 
that  the  boundary  layer  at  SF2  was  driven  pri¬ 
marily  by  SI  induced  by  a  down-front  wind. 

The  structure  of  the  boundary  layer  also  sup¬ 
ports  this  hypothesis.  SI  acts  to  reduce  the  anom¬ 
alously  negative  PV  by  inducing  a  circulation 
that  increases  the  stratification,  thereby  counter¬ 
acting  the  effect  of  the  EBF  (Fig.  4).  Simulated 
boundary  layers  within  symmetrically  unstable 
fronts  are  simultaneously  stratified  and  turbulent 
(5),  in  contrast  to  those  outside  of  fronts,  which 
are  generally  well  mixed.  Indeed,  the  observed 
density  profiles  within  the  front  (Fig.  3B)  lack 
mixed  layers  and  are  instead  stratified  at  all  depths. 
The  Lagrangian  float  trajectories  repeatedly  cross 
this  stratification,  indicating  that  the  boundary 
layer  at  SF2  is  both  turbulent  and  stratified  (30). 

Although  SF1  exhibits  elevated  EBF  and  dis¬ 
sipation,  the  thin  (H  ~  10  m)  boundary  layer 
precludes  estimating  PV  and  the  towed  surveys 
barely  cross  the  front,  making  EBF  errors  large. 


Fig.  4.  Structure  of  the 
symmetrically  unstable  front. 

A  wind  blowing  down  the 
frontal  boundary  between 
warm  and  cold  water  induces 
an  Ekman  transport  perpen¬ 
dicular  to  the  wind  and  to 
the  front  This  carries  heavy 
water  from  the  cold  side  of 
the  front  over  light  water 
from  the  warm  side,  which, 
in  the  presence  of  the  fron¬ 
tal  jet  and  lateral  density 
gradient,  acts  to  reduce  the 
stratification  near  the  surface 
and  makes  the  front  unstable 
to  symmetric  instability.  The 
instability  draws  energy  from 
the  frontal  jet,  leading  to  en¬ 
hanced  turbulence,  and  in¬ 
duces  a  circulation  acting  to 
bring  warm  water  to  the  sur¬ 
face  and  cold  water  to  depth, 
thus  counteracting  the  effect 
of  the  Ekman  transport  and 

keeping  the  near-surface  stably  stratified,  with  warm  water  over  cold  water. 
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An  accurate  evaluation  of  the  hypothesis  is  not 
possible  at  SF1. 

Near-inertial  frequency  waves.  Sections  of 
velocity  and  shear  (Fig.  2,  B  and  D)  show  that  the 
above  frontal  processes  are  associated  with  deeper 
structures  suggestive  of  internal  waves.  In  partic¬ 
ular,  the  depth-time  section  of  shear  (Fig.  2D) 
shows  alternating  diagonal  stripes  of  positive  and 
negative  shear  with  upward  phase  propagation 
and  a  period  close  to  the  local  inertial  period  (i.e., 
half  a  pendulum  day:  0.84  days  at  this  latitude). 
The  north-south  component  of  shear  (not  shown 
in  Fig.  2)  is  in  quadrature  with  the  east-west 
component  such  that  the  velocity  vector  rotates 
clockwise  with  approximately  constant  magni¬ 
tude  as  a  function  of  both  increasing  depth  and 
increasing  time.  This  pattern  is  widely  found  in 
the  ocean  and  interpreted  as  the  signature  of  down¬ 
ward  propagating  near-inertial  frequency  internal 
waves  (31).  Given  the  observed  stratification  and 
estimating  the  vertical  wavelength  and  period  of 
the  waves  to  be  200  m  and  0.78  days,  respec¬ 
tively  (based  on  a  least-squares  fit  on  the  shear 
field  of  the  upper  150  m  and  first  2  days),  theory 
predicts  that  the  waves’  downward  energy  flux  is 
~6  mW/m2,  which  is  similar  to  an  estimate  for  the 
energy  input  to  near-inertial  waves  from  the  winds 
of  9  mW/m2  (32),  but  only  about  6%  of  the  ex¬ 
cess  turbulent  dissipation  at  the  SF2.  These  cal¬ 
culations  suggest  that  the  waves  are  probably 
driven  by  the  winds  and  minimally  contribute  to 
the  energetics  of  the  turbulence  within  the  bound¬ 
ary  layer  at  the  front. 

Surprisingly,  however,  the  strong  near-surface 
shear  of  the  sharpest  front  appears  to  be  part  of 
the  deeper  near-inertial  pattern.  The  boundary- 
layer  depth  (Fig.  3E)  also  appears  to  have  vari¬ 
ability  on  roughly  the  same  time  scale;  that  is,  the 
increased  depth  at  days  138.7  and  139.6.  Thus,  it 
is  possible  that  these  inertial  motions  could  play  a 
role  in  the  rapid  confluence  and  difluence  that 
generate  and  dissipate  the  SF,  as  well  as  in  pro¬ 
ducing  its  negative  PV.  We  further  speculate  that 
the  SI  at  the  front  could  feed  energy  into  the  in¬ 
ertial  waves  and  thus  radiate  energy  into  the  ocean 
interior.  Because  the  lateral  scale  of  the  near- 
inertial  motions  is  probably  much  larger  than  that 
of  the  SF,  their  overall  role  in  the  SF  energetics 
could  be  substantially  larger  than  that  implied  by 
the  small  local  flux  density. 

Implications.  Traditionally,  the  upper-ocean 
boundary  layer  is  thought  to  be  driven  by  the 
atmosphere  through  fluxes  of  heat,  moisture,  and 
momentum  (33,  34).  The  observations  presented 
here  break  from  this  paradigm  by  suggesting  that 
lateral  density  gradients  and  their  geostrophic  cur¬ 
rents  can  also  play  a  role  in  boundary-layer  dy¬ 
namics  by  supplying  energy  to  turbulence  at  the 
expense  of  the  circulation  and  permitting  strat¬ 
ification  and  turbulence  to  coexist.  Therefore,  the 
greatly  enhanced  boundary-layer  turbulence  and 
dissipation  described  here  in  a  very  sharp  Kuroshio 
front  is  likely  an  extreme  example  of  a  process 
that  occurs  much  more  widely  in  the  ocean,  po¬ 
tentially  playing  an  important  role  in  its  dynamics 
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and  energetics.  Furthermore,  these  results  are  con¬ 
sistent  with  recent  theory  on  submesoscale  pro¬ 
cesses  and  thus  encourage  incorporation  of  this 
theory  into  boundary-layer  models.  Such  physics 
is  not  accounted  for  in  present-day  climate  models. 
Fronts  associated  with  the  Kuroshio,  Gulf  Stream, 
and  Antarctic  Circumpolar  Current  are  key  players 
in  the  ocean-atmosphere  climate  system.  Inaccu¬ 
rate  representation  of  the  boundary  layer  and  flow 
energetics  in  frontal  regions  could  thus  substan¬ 
tially  affect  the  predictive  skill  of  climate  models. 
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Activation  of  G  protein-coupled  receptors  upon  agonist  binding  is  a  critical  step  in  the 
signaling  cascade  for  this  family  of  cell  surface  proteins.  We  report  the  crystal  structure  of  the 
A2A  adenosine  receptor  (A2AAR)  bound  to  an  agonist  UK-432097  at  2.7  angstrom  resolution. 
Relative  to  inactive,  antagonist-bound  A2AAR,  the  agonist-bound  structure  displays  an  outward 
tilt  and  rotation  of  the  cytoplasmic  half  of  helix  VI,  a  movement  of  helix  V,  and  an  axial  shift  of 
helix  III,  resembling  the  changes  associated  with  the  active-state  opsin  structure.  Additionally,  a 
seesaw  movement  of  helix  VII  and  a  shift  of  extracellular  loop  3  are  likely  specific  to  A2AAR  and 
its  ligand.  The  results  define  the  molecule  UK-432097  as  a  "conformationally  selective  agonist" 
capable  of  receptor  stabilization  in  a  specific  active-state  configuration. 


G protein-coupled  receptors  (GPCRs)  are 
critical  cellular  signal  transduction  gate¬ 
keepers  for  eukaryotic  organisms.  Some 


GPCRs  are  activated  by  ligands  that  act  as  ago¬ 
nists;  others  are  inactivated  by  inverse  agonists 
and  antagonists.  Efforts  to  elucidate  the  crystal 


structures  of  GPCRs  bound  to  diffusible  ligands 
have  recently  yielded  structures  for  five  class  A 
(rhodopsin-like)  GPCRs:  the  p2  (7-7)  and  fa  (8, 9) 
adrenergic  receptors  (p2AR  and  piAR),  the  A2A 
adenosine  receptor  (A2AAR)  (10),  the  CXCR4 
chemokine  receptor  (77),  and  the  D3  dopamine 
receptor  (12).  All  of  these  structures  display  a 
common  seven-transmembrane  (7TM)  topology 
for  GPCRs,  as  well  as  substantial  variations  in 
functionally  divergent  regions — especially  on  the 
extracellular  side  of  the  receptor,  which  is  re¬ 
sponsible  for  the  recognition  of  a  vast  variety  of 
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